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Abstract. Supernovae are the dominant energy source for driving turbulence within the in- 
terstellar plasma. Until recently, their effects on magnetic field amplification in disk galaxies 
remained a matter of speculation. By means of self-consistent simulations of supernova-driven 
turbulence, we find an exponential amplification of the mean magnetic field on timescales of 
a few hundred million years. The robustness of the observed fast dynamo is checked at differ- 
ent magnetic Reynolds numbers, and we find sustained dynamo action at moderate Rm. This 
indicates that the mechanism might indeed be of relevance for the real ISM. 

Sensing the flow via passive tracer fields, we infer that SNe produce a turbulent a effect which 
is consistent with the predictions of quasilinear theory. To lay a foundation for global mean-field 
models, we aim to explore the scaling of the dynamo tensors with respect to the key parameters 
of our simulations. Here we give a first account on the variation with the supernova rate. 
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We here present new results on our local box simulations of a differentially rotat- 
ing, vertically stratif ied, turbulent interstellar medium threaded by weak magnetic fields 
( Gre ssel et al ] l2008ah . In our model, we apply optically thin radiative cooling and heating 
to account for the heterogeneous, multi-phase nature of the ISM. Improving over existing 
models, we compute a radiatively stable initial solution to avoid the transient collapse 
seen in models applying an isothermal stratification. The central feature of our simu- 
lations is the driving of turbulence via several thousand localised injections of thermal 
energy, which well resemble the kinetics of the supernova feedback. Unlike for artificial 
forcing, the energy and distribution of the SNe are determined by observable parameters. 

In a preceding paper, we have shown that t he turbulence created by SNe does in fact 
exponentially amplify the mean magnetic field ( Gressel et al.ll2008bl ). Here we extend this 
work towards a broader parameter base. The organisation of this article is as follows: in 
Section [TJ we discuss the possible relevance of the discovered effect for realistic Reynolds 
numbers and compare the kin etically driven SN dynamo with the cosmic ray dynamo 
found by lHanasz et al. (|2004lh In Section d we then report on the influence of the 
supernova rate on the measured dynamo parameters. For a short review on mean-field 
modelling of the galactic dynamo, we refer the reader to Elstner et al. (this volume). 



1. Slow versus fast dynamo 

In laminar dynamos, diffusion sets the relevant timescale for magnetic reconnection, 
thus defining an upper limit for the allowable growth rate of the mean magnetic field. 
Because the microscopic diffusivity is usually low, these dynamos are commonly referred 
to as "slow dynamos". Within the ISM, the diffusion time Td = L 2 /r/ (related to the 
microscopic value i] ~ 10 8 cm 2 s _1 ) by far exceeds the Hubble time. This means that 
the field amplification mechanism in galaxies needs to be a "fast dynamo" in the sense 
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Figure 1. Evolution of the regular (a) and fluctuating (b) magnetic field strength for different 
Rm. For clarity, the ordinates of the different models have been offset by an order of magnitude 
each. In panel (c), we compare the growth rates for the turbulent (diamonds) and regular 
(triangles) magnetic field. The unconnected data points to the right correspond to a run with 
rj = and provide an indication for the level of numerical diffusivity. 



that it works on a timescale different than Td ( Lazarian fc Vishniac 19991 ). Accordingly, 
one assumes that the galactic dynamo is determined by some sort of effective turbu- 
lent diffusivity rj t . Although the dynamo may still be limited by topological changes 
via reconnection, this process can be considerably faster due to the much higher value 
of rjt compared to rj. Turbulent reconnection has recently been studied numerically by 
Otmianowska-Mazur, Kowal & Lazarian (this volume). 

According to the definition of Td, in the laminar case one expects the dynamo growth 
rate to increase with rj. The picture, however, changes significantly as soon as the 
Reynolds number is high enough to allow for developed turbulence. Structures of integral 
length scale L are now efficiently broken down to the Kolmogorov microscale where the 
atomic diffusion takes over. Varying the microscopic value of rj does only change the 
extent of the inertial range towards higher wavenumbers but is no longer relevant to the 
large-scale flow. This implies that a turbulent dyna mo should be insensitive to v ariations 
in 77 as soon as a critical value Rm c is exceeded (cf. Ihazarian fc Vishniadll999T ). 

In Fig.[T] we plot the evolution of the mean and turbulent magnetic field (normalised to 
the kinetic energy) for different magnetic Reynolds numbers Rm — L 2 qfl/r/. The values 
for Rm are obtained by varying rj while keeping the rotation rate O fixed. The kinematic 
viscosity v is furthermore adopted to keep the magnetic Prandtl number Pm = v /rj = 2.5 
constant. Irrespective of the value of Rm, we observe a nice and steady exponential 
amplification of both the regular and turbulent fields. To estimate the influence of the 
finite grid resolution, we have performed a fiducial run (dark grey lines in panels (a) 
and (b) of Fig. [I} at double the grid spacing for Rm = 10, 000. The obtained values are 
consistent with the higher resolved run and provide a first indication that the simulation 
results are reasonably converged at this level of dissipation and below. 

In panel (c) of Fig. [TJ we compare the growth rates for the turbulent (diamonds) and 
regular (triangles) magnetic field as a function of the magnetic Reynolds number. The 
unconnected data points to the right correspond to a run with rj = 0, i.e., the case 
of (formally) infinite Rm. As can be seen from a comparison with these points, above 
Rm ~ 10, 000 we are limited by the finite value of the inherent numerical diffusivity of 
our code, i.e., better resolved runs become mandatory to study the regime of higher Rm. 
With respect to the reliably converged runs, we observe growth rates that increase with 
Rm - suggesting that the effect is of genuinely turbulent nat ure. 



In t heir models of the cosmic-ray-driven buoyant instabilitv. lOtmianowska-Mazur et al 



(2007) observe that their dynamo crucially relies on the presence of a "microscopic'' 
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Figure 2. Vertical profiles 
of the mean flow u z (light) 
and the turbulent velocity u' 
(dark colour) at quarter, half, 
and full supernova rate a/ao- 



diffusivity n. The authors actually find the efficiency of the field amplification to scale 
with this parameter. Moreover, they state that no scale separation is manifest in the 
spectra of their simulations. In this respect, it remains disputable whether the high 
value for r\ might rather be interpreted as an effective turbulent dissipation rj t . If so, 
the simulations would have to be regarded as large eddy simulations, i.e., simulations 
assuming turbulent effects rather than having them emerge from first principles. 



2. The dependence on supernova rate 

Quasilinear theory (jKrause fe Radlerlll98fj) is a powerful tool in predicting the dynamo 
a effect from the underlying kinetic structure of the turbulent flow. Unfortunately, as far 
as galaxies are concerned, little has been known with certainty about the vertical profiles 
of the turbulent velocity u'(z) and the mean flow u z (z). In par ticular, u z had to be 
ignored in analytical derivations (see e.g. Frohlich fe Schultzll996t ) assuming hydrostatic 



equilibrium. In that regard, direct simulations can give important new insights, especially 
if one is interested in the dependence on factors like the supernova frequency. 

2.1. Vertical structure & wind 

In the following, we measure the supernova rate a in units of the galactic value <7o = 
30 Myr -1 kpc -2 , representative of type-II SNe. All models include sheared galactic ro- 
tation with \qD,\ = 100 km s -1 kpc -1 . The time averaged turbulent and mean velocity 
profiles of three models with ct/<jq = 0.25, 0.5, and 1.0 are depicted in Fig. 03 The ver- 
tical structure of the turbulent velocity u' shows a distinct M-shape, which peaks at 
±1 kpc. The positive gradient of u' in the central disk strongly suggests an inward trans- 
port of the mean magnetic field. The in ner part of the profiles is similarly shaped as 
the ones obtained from MRI turbulence (Dzi ourkevitch et al . 2004; Piontek fe Ostriken 



20071) but considerably steeper. Crudely extrapolating the fall-off in u'{z), we estimate 



that the MRI might become important in maintaining the observed velocity dispersions 
above galactic heights of \z\ ~ 3 kpc. 

While the overall amplitude of the turbulence increases with the SN intensity, its 
gradient is only weakly affected. This can be seen in Fig. 02 where we plot the scaling of 
the fitted velocity gradients. Unlike the turbulent velocity gradient du'/dz, which goes 
into saturation for a ~ 0.1 ctq, the wind profile shows a distinct scaling with the supernova 
frequency. We thus estimate the wind profile from kinetic feedback as 

0.4 



a 



u z (z) ~ 15-kms- 1 — — , (2.1) 

\aoJ l.kpc 

which, of course, neglects the characteristic modulation of the mean flow within the 
V-shaped region of u'(z), where the kinetic pressure counteracts the thermal pressure. 
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Figure 3. Vertical velocity 
gradients as a function of the 
supernova rate a/ao. The val- 
ues are derived from the inner 
disk for u , and from the full 
domain for u z (cf. Fig. [2} ; the 
overplotted regression shows 
a logarithmic slope of 0.4. 
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Figure 4. Coefficients of the dynamo a and fj tensor as obtained with the test-field method. 

2.2. Dynamo tensors & growth rates 

One main focus of our work is the derivation of mean-field closure parameters from direct 
simulations. We here assume a standard parameterisation of the turbulent EMF 



Si 



a ijBj — Vij £ jkldkBi 



i,j e{R,<f)},k = z, 



(2.2) 



which enters the mean-field induct ion equation via V x £ . To m easure these coefficients, 
we apply the test-field method of ISchrinner et al.l (|2005l I2007T ) . From the thus derived 
profiles a(z) and fj(z), we compute vertically averaged amplitudes. In Fig. [H we plot 
these integral mean values as a function of the supernova rate a. Notably, for each of the 
data points we had to perform a separate 3D simulation, covering a few hundred million 
years to obtain reasonable statistics. Despite the strong scatter in the plotted values, we 
observe a quite robust scaling with respect to a. An exception to this are the coefficients 
which display irregular scaling in some of the runs and will be subject to further 
investigations. 

The main result from Fig.|4]is that all coefficients (except ana) scale in a similar way. 
In particular, this implies that the dynamo number 



(2.3) 



remains approximately constant with a. Furthermore, as we have already demonstrated 
in an earlier work, the effects of the off-diagonal elements ancf> and a^n (which are 
responsible for the diamagnetic pumping) are approximately balanced by the mean flow 
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Figure 5. Growth rate t^ 1 of 
the (vertically integrated) reg- 
ular field (B) (diamonds) and 
the turbulent field ( B' } (trian- 
gles) as a function of a/a . 

u z . Comparing the slope shown in Fig. [3] with the ones in the second panel of Fig 2J we 
see that this balance of forces is roughly independent of a and therefore seems to be of 
rather general nature. 

With the constant ratio of the a effect versus the turbulent diffusion rj t on one hand, 
and the mutual balance of the turbulent pumping and mean flow on the other hand, we 
expect the overall growth rates to be largely insensitive to the applied supernova rate. 
For sufficiently high values of a this is indeed the case, as can be seen from Fig. [5] above 
- we want to point out that these growth rates are obtained from the direct simulations 
and hence do not depend on the a prescription. With constant growth rates over more 
than one order of magnitude in star formation activity, the kinetically driven SN dynamo 
proves to be rather universal. Accordingly, it may well be applicable to a wide range of 
scenarios - covering dwarf galaxies as well as starbursts. 




3. Conclusions 

We have presented recent simulation results on the variation of the key dynamo pa- 
rameters with the supernova frequency. Yet more work is needed to understand how the 
vertical structure of the galactic disk is linked to the emergence of the dynamo effect and 
the associated transport processes. As a primary question, one is of course interested 
in finding a simple explanation for the observed scaling relations - which are somewhat 
flatter than the square-root of the supernova rate a. 

From quasilinear theory, it is understood that the dynamo effect depends on vertical 
gradients both in the density p and the turbulence intensity u' . The disk structure, on 
the other hand, is determined by the balance of the kinetic pressure from the forcing 
and the gravitation of the stellar component - here equipartition arguments can yield 
a first approximation. To connect kinematic properties with dynamo parameters, we 
however need to specify the correlation time t c of the turbulent flow. By comparing direct 
simulations with SOCA predictions, it will eventually become possible to determine how 
this quantity depends on the supernova frequency. 

To connect the measured dynamo parameters a and f] with the growth rates observed 
in the simulations - and thereby check the applicability of the approach - more mean- 
field modelling is necessary. We are, however, confident that the classical framework of 
MF-MHD remains a valuable tool in the quest of understanding galactic magnetism. 
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